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ABSTRACT 

We interpret the observed X-ray morphology of the central part of the 
Crab Nebula (torus + jets) in terms of the standard theory by Kennel and 
Coroniti (1984). The only new element is the inclusion of anisotropy in the 
energy flux from the pulsar in the theory. In the standard theory of relativistic 
winds, the Lorentz factor of the particles in front of the shock that terminates 
the pulsar relativistic wind depends on the polar angle as 7 = 70 + 7mSin^^, 
where 70 ~ 200 and 7^ ~ 4.5 x 10^. The plasma flow in the wind is isotropic. 
After the passage of the pulsar wind through the shock, the flow becomes 
subsonic with a roughly constant (over the plerion volume) pressure P — |ne, 
where n is the plasma particle density and e is the mean particle energy. 
Since e ~ ■ymc^, a low-density region flUed with the most energetic electrons is 
formed near the equator. A bright torus of synchrotron radiation develops here. 
Jet-like regions are formed along the pulsar rotation axis, where the particle 
density is almost four orders of magnitude higher than that in the equatorial 
plane, because the particle energy there is four orders of magnitude lower. The 
energy of these particles is too low to produce detectable synchrotron radiation. 
However, these quasi-jets become comparable in brightness to the torus if 
additional particle acceleration takes place in the plerion. We also present the 
results of our study of the hydrodynamic interaction between an anisotropic 
wind and the interstellar medium. We compare the calculated and observed 
distributions of the volume intensity of X-ray radiation. 

Key words: plasma astrophysics, hydrodynamics and shock waves. 
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INTRODUCTION 

The Crab Nebula is one of the most interesting and best-studied sources in the sky. 
This object was observed over a wide wavelength range: from radio to gamma-rays with 
a photon energy of 50 TeV (Aharonian 2000; Hester 1998; Shklovskii 1968). However, 
advances in observational astronomy provide new data on the Crab Nebula. In the last 
decade, progress in the technology of X-ray telescopes has allowed the Crab Nebula to be 
observed with an angular resolution comparable to the angular resolution of ground-based 
optical telescopes. As would be expected, this gave completely new information on the 
structure of the central part of the nebula. The data obtained on the Chandra X-ray 
observatory arouse particular interest. 

The Chandra observations show that the central part of the nebula consists of 
two components in the soft X-ray emission: a toroidal structure surrounding the pulsar 
PSR 0531 + 21 and two jet-like features located perpendicular to the torus and emerging 
from the pulsar (Weisskopf et al. 2000). Interestingly, such a structure of the plerion 
central region is also observed in the Vela pulsar (Pavlov et al. 2001) and in the supernova 
remnant GO. 9 + 1 (Gaensler 2001). In this paper, we focus our attention on the Crab 
Nebula primarily because of the parameters of the relativistic plasma flow (pulsar wind) 
from this pulsar. Weak magnetation of the wind from PSR 0531 -|- 21 allows the problem of 
its interaction with the interstellar medium to be simplified to an extent that the intensityy 
distribution of synchrotron radiation in the plerion can be easily estimated and compared 
with the observed one. 

The first impression that arises when studying the Chandra images of the Crab Nebula 
is that the toroidal structure surrounding the pulsar suggests the presence of an accretion 
disk around the nebula. However, this interpretation of the observed picture is incorrect 
for obvious reasons. First, there is no independent evidence for the existence of a second 
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companion and an accretion disk around the Crab pulsar. Second, the characteristic size of 
the toroidal structure itself, ~ 10^^ cm, rules out the possibility of the disk interpretation 
of the observed picture. 

The situation with the jet-like features is more comphcated. Jets are observed from 
many Galactic YSO (Livio 1999), SS 433 (Cherepashchuk 1998), superluminal sources 
(Mirabel and Rodrigues 1998) and AGN (Urri and Padovani 1995). A direct analogy 
between the jets from these objects and those observed in the Crab Nebula suggests itself. 
The assumption that the pulsar itself ejects coUimated plasma flows seems reasonable 
enough. However, it is most likely incorrect. 

The Crab Nebula is a typical plerion — a bubble of relativistic particles frozen in a 
magnetic field. This bubble is formed when the flow of supersonic relativistic plasma (pulsar 
wind) ejected by the pulsar interacts with the interstellar medium. Since the wind itself is 
cold, it remains unobservable thus far (see, however, Bogovalov and Aharonian 2000). After 
passing through the wind-terminating shock, the particles are isotropized and begin to emit 
synchrotron photons over a wide electromagnetic spectral range, producing the observed 
plerion radiation; hence the fundamental difference between the jet-like features observed 
in the nebula and the actual jet flows observed in the Universe from other objects. The 
latter are the supersonic coUimated flows ejected from the source. They are characterized 
by termination when interacting with the interstellar medium to produce a shock and the 
so-called lobes (Ferrari et al. 1996). The jet-like features themselves in the Crab Nebula 
are formed behind the shock that terminates the pulsar wind. The plasma flow in them is 
definitely subsonic. Therefore, the physics of this phenomenon undoubtedly differs from the 
physics of the processes that give rise to astrophysical jets. 

Since the observed structures in the Crab Nebula result from the interaction of the 
pulsar wind with the ambient medium, the pattern of this interaction must be studied 
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to understand their nature. The interaction of the wind from the Crab pulsar with the 
interstellar medium has been analyzed by many authors. Rees and Gunn (1974) and Kennel 
and Coroniti (1984) first gave important constraints on the parameters of the pulsar wind 
immediately in front of the shock. The calculated plerion expansion velocity, luminosity, 
and synchrotron radiation spectrum (from optical wavelengths to X-rays) agree with the 
observed ones if the wind consists of electrons and positrons with a Lorentz factor of 
~ 3 X 10^ and if almost all of the pulsar rotational losses are transformed into the particle 
kinetic energy so that the ratio of the electromagnetic energy flux to the particle kinetic 
energy flux is u = 3 x 10~^. For such wind parameters, we can also naturally explain the 
gamma-ray emission from the Crab Nebula above 10 GeV, which is generated by the inverse 
Compton scattering of the same electrons that generate the synchrotron radiation (Atoyan 
and Aharonian 1996; de Jager and Harding 1992). 

The success of the theory by Kennel and Coroniti (1984) was achieved through a 
signiflcant simpliflcation of the problem. This theory assumes the problem to be spherically 
symmetric. As long as the analysis was restricted to the integrated characteristics of the 
radiation from the Crab Nebula (spectra, luminosity), this limitation was not fundamental 
in nature. However, the observed X-ray morphology of the Crab Nebula cannot be explained 
in terms of this theory. Clearly, the Crab Nebula is not spherically symmetric. The more 
realistic pattern of interaction between an anisotropic pulsar wind and the interstellar 
medium must be considered. Here, we made the flrst step in solving this problem. At this 
stage, we do not set the goal of developing a full-blown theory of the interaction between 
an anisotropic, magnetized pulsar wind and a homogeneous interstellar medium. This 
is not yet possible. Here, we determine the pattern of pulsar-wind anisotropy by using 
the results that have been obtained in pulsar physics in recent years. We also perform 
a semiquantitative analysis of the result of the interaction between such a wind and the 
interstellar medium, including an estimation of the synchrotron radiation. In the end, we 
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wish to understand whether the level of anisotropy in pulsar winds that follows from the 
pulsar theory is enough to explain, at least in general terms, the structure of the central 
part of the Crab Nebula observed on the Chandra observatory. 

THE PULSAR WIND FROM PSR 0531 + 21 

The integrated characteristics of the Crab Nebula can be naturally explained in terms 
of the theory by Kennel and Coroniti (1984) if the wind magnetization parameter is 
(J = 3 X 10~^ and the wind Lorenz factor is 7 ~ 3 x 10^ immediately in front of the shock 
(Kennel and Coroniti 1984), although other wind parameters cannot be completely ruled 
out either (Begelman 1998). The Kennel-Coroniti theory naturally accounts for the source 
spectrum over an unprecedentedly wide wavelength range, from optical to hard gamma-ray 
emission (fifteen orders of magnitude in wavelength!). No theory in astrophysics can boast a 
similar success. However, there is one problem in this theory that spoils the overall picture. 
It is not yet clear how the pulsar PSR 0531 + 21 produces the relativistic wind with such 
parameters. This remains one of the key puzzles in the physics of radio pulsars. 

The problem is that all the currently available theories of particle acceleration and 
plasma formation in pulsar magnetospheres [the polar-cap theory (Arons 1983; Daugherty 
and Harding 1996) or the outer-gap theory (Romani 1996; Cheng et al. 2000)] are capable 
of explaining how the dense plasma that produces a relativistic particle wind is formed. 
However, this plasma carries a negligible fraction (~ 10~^) of all rotational losses from the 
Crab pulsar. The entire energy flux from the pulsar is concentrated in the electromagnetic 
field carried away by the wind, which corresponds to the wind magnetization parameter 
a ~ 10^. The wind itself has a modest Lorentz factor ~ 200 near the radio- pulsar light 
cylinder (Daugherty and Harding 1996). It is yet tbe clarified through which processes 
almost the entire electromagnetic energy fiux is transformed into the wind-particle kinetic 
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energy on the way from the hght cyhnder to the shock front, although substantial efforts 
were spared to solve this problem (Coroniti 1990; Lyubarsky and Kirk 2001). 

There is no need to know the wind acceleration mechanism to determine the 
energy-flmux distribution in the wind. The fact that the conservation of the energy flux in 
the wind holds in any case is suffice. The point is that the electromagnetic energy flux in 
the winds from radio pulsars propagates along streamlines. The plasma kinetic energy flux 
also propagates along these streamlines. This conclusion is based on MHD models of the 
winds from axisymmetrically rotating objects (Okamoto 1978). However, it has recently 
been shown to be also valid for obliquely rotating objects where the flow beyond the light 
cyhnder is concerned (Bogovalov 1999). This ensures that the total energy flux per particle 
is conserved along a given streamline. The kinetic energy flux per particle, in units of mc^, 
is 7. The electromagnetic energy flux per particle, in units of mc^, is 



s ^ — [E X B]/{nvmc^'j), (1) 

where E is the electric held, B is the magnetic held, n is the particle density in the intrinsic 
frame of reference, and v is the wind velocity. The sum 7 + s depends on the streamline but 
is conserved along it. It follows from the solution of the problem on the structure of the 
wind from an oblique rotator (Bogovalov 1999) that the energy flux in the wind sufficiently 
far from the light cylinder may be considered to be azimuthally symmetric, although the 
electromagnetic field itself is not azimuthally symmetric in this case (see Bogovalov (1999) 
for details). In this notation, the magnetization parameter is cr = s/7. The conservation of 
the total energy fiux ensures that 



7o + So = 7 + s- (2) 
The subscript '0' marks the values near the hght cylinder. Since s <^ 7 in front of the shock 
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that terminates the pulsar wind, the Lorentz factor of the preshock plasma may be assumed 
to be 7 = 7o + sq. Thus, the dependence of the particle energy on the streamline along 
which the particles move is determined by their initial energy and the initial distribution 
of the electromagnetic energy flux. We make the only assumption. Assume that the 
initial distribution of the electromagnetic energy flux does not depend on (or is almost 
independent of) the wind acceleration. This assumption holds true in all cases if the wind is 
accelerated sufficiently far from the light cylinder in the supersonic flow region. Then, the 
field distribution in the magnetosphere and, hence, the initial electromagnetic energy fiux 
do not depend on what happens downstream of the magnetosonic surface. This is because 
no MHD signal can penetrate from the supersonic fiow region into the subsonic fiow region 
and affect the fiow in this region. If this is the case, then it will suffice to determine Sq, 
provided that there is no acceleration. 

Numerical and analytical calculations of the relativistic plasma fiow show that the 
plasma magnetic coUimation is neghgible for Lorentz factors 7 > 200 (Beskin 1998; 
Bogovalov and Tsinganos 1999; Bogovalov 2001). The pulsar wind may be assumed 
to spread out radially. Since the relation E — ^^^^J^ Bp (Mestel 1968), where r is the 
distance to the pulsar, 9 is the polar angle, is the pulsar angular velocity, and Bp is the 
poloidal magnetic field in the wind, holds between the electric field in the wind and the 
poloidal magnetic field, the condition for the magnetic field being frozen in the plasma 
E+l[Ex B]^0 takes the form 

r sin 9D,Bp + B^Vp — v^Bp. (3) 

At r ^ c/fl, the plasma angular velocity tends to zero, because the angular momentum 
of the plasma particles (rsin(^)v<^) is limited above. We then derive a simple expression 
for the toroidal field in the wind far from the pulsar, B^ — '•^"1^^ r Therefore, the 
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electromagnetic energy flux per particle (in units mc^) is 



Vp J 47rnomc^7o 



So = sin^^ — (4) 



We see that when the relativistic wind spreads out radially and uniformly, the Lorentz 
factor 7i of the preshock wind particles must have a latitude dependence of the form 

7i = 7o + 7m sin^ 9, (5) 



where 7™ = . % ■ This is the maximum Lorentz factor of the preshock wind 

particles. To be consistent with the theory by Kennel and Coroniti (1984), it must be of 
the order of 3 x 10^. 

Note that expression (5) for the particle Lorentz factor immediately follows from the 
MHD theory of magnetized winds from rotating objects and is almost model-independent. 
It only assumes that the particle flux from the pulsar is isotropic. Clearly, this assumption 
does not severely restrict the range of applicability of our results. For the standard 
parameters 70 = 200 and 7,^ ~ 3 x 10^, the Lorentz factor changes with latitude by four 
orders of magnitude. Even if the particle flux changes with latitude by several times (or 
several tens of times), this does not change the overall dependence. Anyway, the most 
energetic particles will be near the equator and their Lorentz factor will be higher than the 
Lorentz factor of the particles at the rotation axis by several orders of magnitude. 

Below, for our calculations, we assume the wind flow to be radial with an isotropic 
mass flux. The plasma density in the intrinsic frame of reference in such a wind is 



_ N 



Here, is the rate of particle injection into the nebula. The toroidal magnetic fleld has the 
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distribution 

Boro sin 9 

B^ = , (7) 

The quantity Bq can be determined from the condition a — 3 x 10~^. We ignore the 
poloidal magnetic field, because it is much weaker than the toroidal magnetic field ahead of 
the shock front. The plasma Lorentz factor is given by expression (5). Our objective is to 
determine (at least qualitatively) how a wind with the above parameters interacts with a 
homogeneous ambient interstellar medium. 



THE APPROXIMATIONS 

The problem of the interaction between a highly anisotropic, magnetized relativistic 
wind and a homogeneous ambient medium has no analytic solution. Even obtaining a 
numerical solution seems problematic so far. Below, wc make two reasonable simplifications 
that will allow us to answer the questions of interest by using simple mathematics. 

(1) The approximation of a hydrodynamic interaction. In the special case of the 
Crab Nebula, as a first approximation, we may disregard the magnetic-field effect on the 
postshock plasma dynamics. This approximation seems reasonable, because the pulsar 
wind from PSR 0531 + 21 is weakly magnetized. Recall that the ratio of the preshock 
Poynting fiux to the plasma kinetic energy fiux is cr = 3 x 10~^. Although the magnetic field 
increases in strength by a factor of 3 after the shock passage, the ratio of magnetic pressure 
to plasma pressure is small up to distances approximately equal to five shock radii (see 
Kennel and Coroniti 1984). Further out, the magnetic pressure is higher than the plasma 
pressure and it cannot be ignored. However, the region of Srgh completely suits us. It is in 
this region that the most interesting features of the central part of the Crab Nebula are 
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formed: the X-ray torus and the jet-hke features. Therefore, below, the wind interaction is 
considered as a purely hydro dynamic one. We will determine the magnetic-field evolution 
from the induction equation with a given distribution of the plasma velocity field. Note 
that although the magnetic-field effect on the plasma dynamics immediately behind the 
shock wave is marginal, Lyubarsky (2002) attempts to explain the observed jets in the 
Crab Nebula as resulting from a magnetic compression of the wind after the shock wave. 
Below, we show that these jets are formed in the hydrodynamic approximation without any 
involvement of the magnetic field. 

(2) The quasi- stationary approximation. Strictly speaking, the interaction of a 
supersonic wind from the central source with the interstellar medium is not stationary. It 
is schematically shown in Fig. 1. The source continuously injects new particles into the 
nebula and the nebula size increases with time. However, in a bounded region of space 
near the shock, the flow may be considered to be steady, provided that the nebula size 
is much larger than the characteristic size of the shock. This is easy to understand from 
simple considerations. Assume that the source injects particles with characteristic energy e 
into the nebula at a rate N . The nebula size Rn at constant external pressure Pext is then 
determined by the relation 7Vie/(47r/3i?^) = Pext, where t is the source operation time. We 
see from this relation that the nebula radius increases as t^/^. This means that the velocity 
of the nebula outer rim decreases with time; therefore, the entire flow may be considered 
in the limit i — > oo as steady with the boundary condition at infinity ^ for r ^ oo 
and P — Pext- Thus, the condition for applicability of the quasi-stationary approximation 
is Rn ^ ^sh, where r^-^ is the characteristic radius of the shock front. For the Crab Nebula, 
i?„ w 2 pc and the shock radius is Tgh = 0.1 pc (Kennel and Coroniti 1984). Consequently, 
the postshock flow for the Crab Nebula within Srgh may actually be considered to be steady. 
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THE INTERACTION OF A HIGHLY ANISOTROPIC WIND WITH THE 

INTERSTELLAR MEDIUM 

Basic Simplifications 

The problem of the interaction between a supersonic, highly anisotropic plasma flow 
and a homogeneous medium has no exact solution so far. Therefore, below, to estimate 
the observed effects that must arise during such an interaction, we proceed as follows: 
flrst, we use a highly simplified interaction model to calculate the volume luminosity of 
the plerion produced by a Crab-type pulsar; subsequently, we qualitatively consider how 
our simplifications affect the results of our calculations by using, in particular, the results 
obtained for weak anisotropy in Appendix A. 

We use the following simplifications to estimate the volume luminosity of the Crab 
Nebula: 

(1) Since the postshock plasma fiow is subsonic, with the plasma velocity tending to 
zero when moving downstream of the shock, the plasma density along a streamline may 
be assumed to be constant (Landau and Lifshitz 1986). The subsonic motion also imphes 
an approximate equality of the pressure in the plerion. Previously, Begelman (1992) used 
this approximation to describe the postshock flow. Let us estimate the accuracy with 
which these conditions may be considered to be satisfled in our speciflc case. The pressure 
variation in the plerion is AP/P ~ Wu'^/P. Here, W = AP and u ~ l/\/8, because 
the plasma velocity is v — \c immediately behind the shock and then rapidly decreases. 
Therefore, in the worst case, we have AP/P ~ 0.5. This error completely suits us, because 
below, we are concerned with the variations in mean particle energy e and plasma density n 
across a streamline, which are four orders of magnitude; the pressure is related to these 
quantities by P = |ne. Against the background of such variations, the pressure variations 
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of 50% are of no fundamental importance. 

(2) We assume the postshock streamlines to remain radial without bending at the 
shock and use the conditions for a perpendicular shock to determine the postshock plasma 
parameters. 



The Shape of the Shock Front 

To determine the shape of the shock front, we use the shock-adiabat relation (A. 18) for 
an oblique shock wave. After the passage of the shock front, the plasma on each streamline 
adiabatically decelerates; far from the shock front, its velocity tends to zero and the plasma 
pressure comes into equilibrium with the external pressure Pext- The relationship between 
the plasma parameters immediately after the shock and Pext is given by the Bernoulli 
equation for a relativistic plasma 



4^2— = 4^ = 7imcl (8) 

Here, we use the constancy of the Bernoulli integral at the shock front and the relativistic- 
plasma approximation e = 3p, which holds good behind the shock. Below, the subscripts '1' 
and '2' denote the preshock and postshock quantities, respectively. 

It is easy to find from these relations and from the adiabatic flow condition (A. 5) that 



ext ^ I ^ \ 



3 / \ <5-i 2-S 

7i TT-i mc^ 2 



it- (9) 



The location of the shock front could be determined from this equation if 7|| were known. 
Since the problem is azimuthally symmetric, this quantity is known only at points on the 
equatorial plane and on the rotation axis. The shock front crosses them at a right angle. 
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Below, we consider a simplified case by assuming, for simplicity, that 7|| = 1 everywhere. 



We then obtain for the shock radius 



3 / 3 ^imc'^N 



(10) 



^2 1,^8^ 47rPextc' 



In the limit 7^ 3> 70 of interest, the shock radius may be assumed to be Vgh = ^cql sin 9 



everywhere, except for a narrow interval of angles 9 < yjo/jm- This implies that the shock 
front is generally a torus whose cross section is a couple of contacting circumferences with 
the centers on the equatorial plane in the middle of the distance r^q from the source to the 
shock. 

To estimate an error in using everywhere the approximation 7|| = 1, we must solve the 
problem more accurately. Our analysis, which is beyond the scope of this paper, shows 
that the shock front lies at slightly larger distances than (10). In addition, a system of two 
shocks is formed near the axis for a highly anisotropic wind. Nevertheless, expression (10) 
gives an error in the shock radius within 15%. In this paper, such an accuracy is admissible. 

The Formation of a Toroidal Structure 

As follows from the condition of a constant density along the postshock streamline, the 
plasma velocity Vp on each streamline behaves as 



where rsh{9) is the distance from the pulsar to the shock location at a given angle 9. This 
dependence on r follows from the conservation of mass flux for radial plasma motion. 

On a given streamline, the dependence of the magnetic field on r and 9 follows from 




(11) 
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the frozen-in condition. According to this condition, B/nr — const on a streamline (Landau 
and Lifshitz 1982). For the preshock magnetic field, we use the condition 

(Bo^smdf -y^sin^^ 
Anriimc^'-fi 71 

It implies that the ratio of the Poynting fiux density to the plasma kinetic energy fiux 
density is everywhere equal to the same value a, except for a narrow region near the 

rotation axis where the toroidal field must vanish. We inserted the additional factor ^ 

7 

in the right-hand part of Eq. (12) to take into account this circumstance. It follows from 
expression (9) that the preshock magnetic field is 



Bi — Bo—sm9 — J 4(T7r7^7inimc^ sin 9. 



(13) 



Given that the magnetic field increases in strength by a factor of 3 after the passage of a 
strong shock, we obtain 



s 

iS-1 



B-i = 3. — Pext sm^ — . 

\ 7 3 V 3 y rsh 



(14) 



At angles 9 ^ vlohm-, the expression for the magnetic field takes the form 



B2 = Seq — , (15) 

where is the postshock equatorial magnetic field. We see that the field immediately 
behind the shock is everywhere the same, except for a narrow region near the rotation axis. 
Outside this region, the field more rapidly increases with distance from the pulsar at high 
latitudes. The field linearly increases until the magnetic energy density becomes equal to 
the plasma energy density. Subsequently, the field begins to decrease with increasing r 



-16- 



(Kennel and Coroniti 1984). The fact that the hnear dependence extends to Srgh (Kennel 
and Coroniti 1984), within which the toroidal structure is formed, will suffice. 

To calculate the synchrotron radiation, we assume, as in Kennel and Coroniti (1984), 
that the following power-law particle spectrum is formed behind the shock: 

F{i, r = Tsh) = A J 77(7 - 7min)r?(7m£« - 7), (16) 

\ 7min / 

where A is the normalization factor determined from the particle injection rate on a given 
streamline, ^min is the minimum particle energy in the spectrum, ^ymax is the maximum 
particle energy in the spectrum, and ri{x) is the step function equal to unity at a; > 1 and 
zero at X < 1. 

The mean energy of the chaotic particle motion in this spectrum must correspond to 
the mean postshock particle energy determined from the shock adiabat. It follows from this 
condition that 



^ TiW (c^ 

72 (a- 1) ^ 

Here, we use the fact that although the cutoff energy of the spectrum, 7inax ~ 5 x 10^ 
(Atoyan and Aharonian 1996), is finite, it is much larger than 7inin < 5 x 10^. 

The plasma moves behind the shock as a whole at velocity (11). The evolution of the 
particle spectrum during this motion is described by the transport equation 

dF d 

The rate of change in the Lorentz factor of the chaotic particle motion, 7, is generallyly 
determined by synchrotron losses and plasma heating thrgh adiabatic compression during 
the wind deceleration. In our case, the adiabatic changes in particle energy can reach 24% 
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of the particle energy. However, we disregard these changes here, because this is a clear 
excess of the accuracy under our assumptions about the postchosk plasma dynamics. 

The solution of Eq. (18), provided that the distribution function at the shock matches 
function (16), is 



Tmin \ T) 



mm , 



In this expression, the function 

describes the degradation of the particle energy through synchrotron losses. The 
function n{r, 9) is the emitting-electron density in the observer's frame of reference. The 
electron density in the plerion on the streamline with 9 follows from expression (6): 



^(^^) = 7rZT77iT- (21) 



3iV 
47r^(^) ' 

Here, we took into account the fact that the density at the shock increases by a factor of 3. 
The synchrotron flux is determined by the convolution of spectrum (19) with the spectral 
distribution of the radiation from an individual electron (Landau and Lifshitz 1973). The 
input parameters for our calculations are the particle injection rate A'" = 5 x 10^^ part. s~^, 
the maximum Lorentz factor of the wind at the equator 7^ = 4.5 x 10^, and the shock 
location at the equator; we correlated the latter with the inner ring of the toroidal structure, 
which is located at = 4.3 x 10^'' cm corresponding to a distance of 14" (Weisskopf et 
al. 2000). Actually, this implies that the external pressure Pext was chosen so that the shock 
at the equator was located on the inner ring of the toroidal structure. 
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The results of our calculations are presented in Fig. 2. This figure shows the plerion 
volume luminosity at photon energy 400 eV, the characteristic energy at which the Chandra 
observations are carried out (Weisskopf et al. 2000). No integration along the line of sight 
was performed. All geometric sizes were reduced to the pulsar distance, 2 kpc. Here, only 
the plerion cross section in the poloidal plane is shown. We see that the spatial distribution 
of the volume luminosity has the shape of a torus with the characteristic sizes corresponding 
to the observed ones. The radiation reaches the highest intensity at a distance of ~ 40", 
where the second outer ring in the Crab Nebula is located. The intensity of the synchrotron 
radiation behind the shock is known to gradually rise (Kennel and Coroniti 1984). This rise 
results from a linear increase in the magnetic field behind the shock, which takes place only 
if the preshock wind was weakly magnetized. The further rapid decline in the radiation at 
distances > 60" is attributable to fast synchrotron electron coohng in the growing magnetic 
field. We see that although our calculations are definitely incorrect for r > 5i?sh, this is of 
no importance. When the electrons reach this region, they have already cooled down to an 
extent that they do not produce detectable radiation in the Chandra spectral range. 

In our calculations, we failed to obtain the bright inner toroidal ring. It may well 
be that this brightening cannot be explained in terms of magnetic hydrodynamics in 
principle and that it is attributable to the structure of a coUisionless shock, which must be 
investigates separately (see Gallant 1992). 

Across the equatorial plane, the agreement with the observed distribution is slightly 
worse. The distribution of the calculated brightness across the equator is broader than that 
of the observed one. It may well be that distributions (5)-(7) do not faithfully describe 
the preshock wind parameters. The energy and plasma fluxes in the wind may be more 
concentrated toward the equator than follows from the theory. In our view, however, it is 
too early to draw this conclusion. The observed disagreement may be entirely attributable 
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to our approximations. We assumed the shock to be perpendicular to the streamhne 
everywhere. This is not true outside the equatorial plane, where the shock inclination to the 
streamline decreases while 7|| increases. As the angle of incidence of the shock decreases, 
the mean energy of the chaotic particle motion after the shock passage also decreases, 
as is clearly seen from the relation for the shock adiabat of an obhque relativistic shock. 
Therefore, the volume intensity must decline with height above the equator faster than in 
our case. 

We disregarded the fact that the streamlines must bend toward the equator at the 
shock, as is the case for weak anisotropy under consideration. Allowing for this circumstance 
will cause an increase in the brightness near the equator and its faster decrease under it. 
Finally, we assumed the streamlines after the shock passage to remain straight. In fact, as 
our calculations for low anisotropy show, the streamlines continue to approach the equator, 
at least in some preshock region, even after the shock passage. Thus, all our most important 
approximations used to estimate the volume luminosity of the plerion produced by the Crab 
pulsar lead to the same effect: a broader distribution of the radiation intensity than must 
be in the case of a more accurate calculation. Therefore, so far we have no reason to believe 
that distributions (5)-(7) inaccurately describe the preshock parameters of the wind from 
the Crab pulsar. 



Jet-like Features 

In the model that is a simple generalization of the model by Kennel and Coroniti (1984), 
we failed to obtain something similar to the bright jet-like features observed in the Crab 
Nebula. However, noteworthy is one circumstance that has a direct bearing on the observed 
jets. As we already pointed out above, the plasma pressure in the plerion is roughly the 
same, because the flow is subsonic. The relativistic-plasma pressure is given by the relation 
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p — |n(r, 9)e, where e is the mean particle energy. Since it is proportional to the Lorentz 
factor of the preshock plasma 71, the plasma density in the plerion is 

n(r,^)=ne,^^^, (22) 
7o + 7mSm 9 

where is the equatorial plasma density. Density (22) appears in expression (19) for 
the electron spectrum. We see from this expression that riaxis/^cq = I'm/lo- In standard 
models, 70 ~ 200 and 7^ a; 3 x 10^, implying that the plasma density in the plerion near 
the rotation axis is approximately a factor of 15 000 higher than the equatorial plasma 
density. It is important to note that this fact is a direct result of the energy distribution in 
the wind (5) and depends weakly on any other factors. Thus, if the luminosity were simply 
proportional to the plasma density, then we would observe only an extremely bright jet 
from the pulsar in the plerion with a barely visible torus against its background. However, 
this is not the case. The increase in density toward the rotation axis causes no brightening, 
because this increase is accompanied by the simultaneous decrease in particle energy. As a 
result, the volume intensity decreases toward the rotation axis. 

Thus, the formation of dense but relatively cold jet-like features in the plerion along 
the rotation axis necessarily follows from the theory of the interaction between the Crab 
pulsar anisotropic wind and the interstellar medium. There is only one problem: through 
which processes these features can become bright enough to be observable. A possible 
solution could be the assumption that some particle acceleration takes place not only at 
the shock but also in the entire volume of the plerion. This assumption seems reasonable 
(Begelman 1998) and is supported by the detection of gamma-rays with energy above 
50 TeV from the Crab Nebula (Tanimori et al. 1998). If we assume that an additional weak 
particle acceleration takes place in the entire volume of the nebula, the spectral shape of 
the accelerated particles is the same everywhere, and the number of accelerated particles 



- 21 - 



is proportional to the plasma density in the plerion, then a second radiation component 
proportional to the particle density in the nebula emerges. We see from Fig. 2 that in this 
case, the second radiation component manifests itself in the form of bright jets, with the 
fraction of the accelerated particles being 2 x 10~^ of their local density. 

CONCLUSIONS 

Based on the standard theory by Kennel and Coroniti (1984), we have shown that the 
principal features of the observed X-ray structure in the Crab Nebula can be naturally 
explained by taking into account anisotropy of the energy flux in the winds from radio 
pulsars. For our calculations, we used the fact that the energy flux density in the pulsar 
wind is proportional to sin^ 6. This dependence follows from the expression for the Poynting 
flux Py for axisymmctric winds, Py ~ r"^ siv? OB^, provided that the poloidal magnetic 
field Bp pulled by the wind from the source is isotropic. Such a dependence remains valid 
even for an obliquely rotating source with a uniform magnetic field, although, in general, 
the flow is not steady and axisymmetric in this case (Bogovalov 1999). For a nonuniform 
(in 9) poloidal magnetic field, the 9 dependence of the energy flux density can change. 
However, it is easy to understand that the energy flux density is at a minimum on the 
rotation axis and reaches a maximum at the equator, irrespective of the rotation angle 
and at any reasonable distribution of the poloidal magnetic field in 9. This is because the 
Poynting flux density is proportional to the toroidal magnetic fleld B^p. However, Bip — 
always on the axis. Therefore, = on the rotation axis and Py can only increase when 
moving away from the axis. 

The breakdown of azimuthal uniformity of the poloidal magnetic field leads to the 
additional generation of magnetosonic waves (but not magnetodipole radiation) in the wind. 
However, the energy flux in them is also proportional to the combination r^sin^6'(5p) (at 
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least for a small wave amplitude; Bogovalov 2001a). Therefore, irrespective of the rotation 
angle, Py always reaches a maximum at the equator, except for the exotic case where Bp 
decreases toward the equator faster than sin~^ 9. In other words, the concentration of the 
energy flux density toward the equator is apparently a common property of the pulsar 
winds. 

The formation of a bright X-ray torus is a direct result of this common property of 
the pulsar winds. The interaction between a wind with such anisotropy and the interstellar 
medium also inevitably gives rise to cold subsonic (jet-like) flows in the plerion along the 
rotation axis whose density is almost four orders of magnitude higher than the equatorial 
density. In the standard theory, these flows are invisible in X-rays, because the particle 
energy is too low to produce detectable synchrotron radiation. However, if we assume the 
additional acceleration of a mere 10~® fraction of the particles at each point of the plerion, 
then the radiation from the jet-like features becomes comparable in intensity to the the 
torus radiation and the overall morphology of the plerion becomes similar to that observed 
on the Chandra observatory. It would be natural to assume that similar features detected 
around the Vela pulsar (Pavlov et al. 2001) and in the supernova remnant GO. 9 -|- 01 
(Gaensler 2001) can be explained in a similar way, because the anisotropy in the pulsar 
wind of the type discussed here must be formed in all pulsars. 
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The actual pulsar wind has a strong latitudinal dependence of the particle energy. 
The ratio 7^/70 is of the order of 10^. However, for a qualitative understanding of the 
interaction between an anisotropic wind and the ambient medium, it is of interest to 
consider this interaction for a weak anisotropy. This problem is valuable in that its solution 
can be obtained analytically. 

As was already pointed out above, the approximation of a hydrodynamic interaction 
is invoked to describe the plasma flow. In addition, the problem is axisymmetric. Let us 
introduce the stream function ip. It is related to the physical quantities by 



A. 



A WEAKLY ANISOTROPIC WIND 



nu0 — 



nUr = 



1 dip 
r sin 9 rdO ' 
1 dij 
r sin 9 dr 



(A2) 



(Al) 



This manuscript was prepared with the AAS LM^^X macros v4.0. 
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where n is the intrinsic particle number density and ug are the corresponding components 
of the four- velocity. It is convenient to use the equations for the stream function in spherical 
coordinates (Beskin 1997) 



2UrUg 



re 



u^nur COS 9 + nu0 sin 9 {2ul + ul) — (A3) 



(nrsin^)^ 



r, r,d\nA S 



1, 9 o.rfln^ 



where ul — and 5 = 4/3 is the polytropic index of the relativistic plasma, 



P 



(A4) 



Since the motion is adiabatic, 5" is conserved along a streamline and depends only on ■0. 
The ratio of the Bernoulli integral and S is 



^ " ^ Jio ^^^^ 

It is also conserved along a streamline. Here, w is the thermal function per particle in the 
intrinsic coordinate system and p is the intrinsic pressure. In the ultrarelativistic case, 
w and p are related by the equation of state w — gzjP- For weak anisotropy, we solve 
Eq. (A. 3) using the perturbation theory by assuming that the Lorentz factor of the wind 
depends on the polar angle as 



7 = 7^(1 + Q;sin2e), (A6) 
where a is a small dimensionless parameter and is the Lorentz factor of the wind. 

In the zero approximation, the flow is spherically symmetric, ug = 0, and the following 
identities hold: ^ — and = 0. The right-hand side of Eq. (A. 3) becomes zero and it 
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takes the form 



^lipoee ~ ulnur cos ^ = 0. 



(A7) 



The solution of Eq. (A. 7) is 

i/jo = nsh-Ushr-sh^oll - cos 6*), (A8) 

where rsh,o is the radius of the spherical shock wave and nsh,Ush are the density and the 
radial component of the four- vector immediately behind the shock. According to Landau 
and Lifshitz (1986), Ush — To solve the equation in the first order of the perturbation 
theory, we represent the stream function as 

i/j = nshUshr^hfl (1 - cos 61 + af{r, 9)) , (A9) 
where /(r, 9) is the correction to the first approximation. The equation for it is 



{rishrihUsh) 2 



^2 "sJVV ^ ^ Jr 



(AlO) 



dr] 



+ 



dr] 



where r] — cos9 and Uq, no, 7 are the solution behind the shock in the zero approximation. 

To determine A and we use the shock-adiabat relations for an oblique relativistic 
shock wave. These are derived from Landau and Lifshitz (1986) using the Lorentz 
transformations 



(P2 -Pl)(g2 +7^1) 

(62 - ei)(ei +P2) 



.1/2 



(Aii; 
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l-vl 



{P2 - Pl){('l +P2) 



1/2 



(A12) 



(62 - ei)(e2 

where p is the pressure and e is the internal energy. The quantities with the subscripts 1 
and 2 describe the preshock and postshock states of the wind, respectively, and v± are 
the tangential and normal velocity components, respectively. 

Since the thermodynamic state of the matter in the preshock region corresponds to a 
zero temperature. 



Pi — 0,ei — m(?ni. (A13) 
According to (21, A. 6), the density can be represented as 



m = (1 -asin2^)^,n^ = (A14) 
The function rsh{0) can be represented as a series of Legendre polynomials P^, 



rsh(^) = r,h,o (^1 + tt ^ RmPmicos 6)^ , (A15) 
In front of the shock, we then have 



n. 



^ll-asm'^e-2aJ2RmPm{cos9)], (A16) 



mc^n. 



ei = [l-a sin^ e-2aJ2 RmPm{cos 9) . (A17) 



^sh,0 



In the ultrarelativistic case, the following relations hold: 62 ^ ei and P2 3> ci. They 
allow (A. 11, A. 12) to be simplified. For 5 = 4/3, we have 



-29- 



P2 = ^^ei, ^i = l^l (A18) 

where 7|| = 



The postshock plasma density can be calculated by using the Bernoulli integral 



72 S P2 . 

n2^ni — - — - — . (A19) 
7i 5 - 1 ei 

For weak anisotropy, v\\ <^ 1. We see from the geometry shown in Fig. 3 that 



4— 5"'-^^— Ksfir- 

Since Tgjj ~ a, the first corrections to 7|| begin with the terms proportional to a^. In the 
first order of smallness in a, we may assume that 7|| = 1. The change in the postshock 
thermodynamic quantities when the flow isotropy breaks down results only from a change 
in the spatial location of the shock front. To a flrst approximation, the change in the 
inclination of the shock surface is negligible. 

In that case, 



2 Tlw 2 2 
P2 = n^—^W^C 



3 ^sh,0 



1 + a sin^ 9-2aY, RmPm{cos d) J , (A21) 

ll^l. (A22) 
n2 = 37on^ ^1 - 2q; ^ RmPm(cos 6)^ . (A23) 



Hence, we have for A and S 



S ^ So(l + asm'^e+ + 2/3a^i?„P„(cos^)j , (A24) 
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where and 5*0 are the constants that are of no interest. 



Substituting the expressions for A and 5" in relation (A. 10) yields 



f -2^f 



2r] 



r 



We seek a solution to Eq. (A. 26) in the form 



(A25) 



(A26) 



n 

where Qm{v) the eigenf unctions of the operator 



(A27) 



(l-n')Ql(n)^-m{m + l)QM; 



These are 



(A28) 



QUv) = (1 - VlPM 



(A29) 



Substituting (A.27) in (A.26) yields 



+ 2/3 I - ' 



(A30) 



where x = r/rsh,o- 



The function uo was numerically determined from the equations for a polytropic, 
spherically symmetric postshock flow. This function depends only on one variable x. Since 
the first boundary condition at the shock front is the continuity of ip, 



= 0. (A31) 

The second boundary condition follows from the shock-adiabat relations for an obhque 
shock wave. Behind the shock, 



?"sh",sli,0"sh<9/ 

n2U2e = -a r-^— (A32) 

r sm B or 

We see from Fig. 3 that at </? <S 1 and /3 <S 1, 



— = <^= — ^^TT' V + P^—- (A33) 

U2r sm Ox ox V2± 

In the first order of smallness in a, we have 



It follows from (A. 33, A. 34) at a; = 1 that 



and the second boundary condition for takes the form 

^— 2iC. (A36) 

The general solution of Eq. (A. 30) with the boundary conditions (A. 31, A. 36) is 
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f^'^^^S^^.ft^ + R^ftK (A37) 

where fi"^^ is the solution of the inhomogeneous equation (.30) with the right-hand side 
_i_2 Us-uo ^j^g boundary conditions fi"^\l) = 0, \x=i= 0. /j™'' is the solution 

of the same equation with the right-hand side ;;^!^| (w^7^ — ""^"" ^ and the boundary- 
conditions fi"^\l) — 0, = 2. Since the functions and /s'"^ increase with distance x 
faster than x^, Rm — when x^ oo under any boundary condition. Therefore, for all 
m 7^ and for m 7^ 2, we may assume that Rm = 0. This condition physically means that 
all obstacles that can be in the flow lie far enough from the shock for their influence on the 
flow near the shock to be ignored. 

Figure 4 shows a family of solutions (A. 37) for various R2. All these solutions satisfy 
the boundary conditions at the shock. To single out the only solution, we must formulate 
the boundary conditions far from the shock. It is easy to verify that if there are obstacles 
in the postshock flow, then the postshock perturbation grows faster than ,x^. In this case, 
the condition that there are no additional obstacles behind the shock in the presence 
of anisotropy is the requirement that the perturbation should grow with x more slowly 
than x^, for example, as This will take place if the right-hand side of Eq. (A. 30) tends 
to zero as x increases. In turn, it is easy to find that the right-hand side of this equation 
tends to zero ior R2 — The solution for this value of R2 is indicated in Fig. 4 by the 
solid line. It is the separatrix that separates the solutions tending to -|-oo from the solutions 
tending to — cxd for a; — > 00. 

Although the expansion of / includes the terms with m > 1, we must also determine Rq, 
because this term appears in the expansion of in terms of Legendre polynomials. This 
coefficient can be easily determined from the condition for the shock location being constant 
on the rotation axis. The shock front then takes the form 
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rsh = r,h,o (l + a (^^^o - ^^^2)) = r,h,o (l + | sin^ 6^ . (A38) 

In this solution, the postshock pressure p2 is constant along the front and the density varies 
as 

?T'2 = 37^n^(l — CKsin^ (A39) 

Accordingly, the mean energy of the chaotic particle motion behind the shock front 
62 = Sp2/n2 is 

2 

^2 = -7«,mc^(l + asin^ 9). (A40) 
3 

Thus, as the energy flux density increases along the equator, the shock front is extended 
along the equatorial plane; the pressure along the front is constant, to a first approximation. 
The postshock particle density at the equator is lower and their mean thermal energy is 
higher than those for the particles along the rotation axis. Figure 5 shows streamlines in the 
flow for a — 0.25. We see that, in addition to this, the increase in the energy flux density 
along the equator also causes the bending of the streamlines toward the equator. Initially, 
this takes place at the shock, where, according to the conditions for an oblique shock wave, 
the streamlines are pressed to the shock. Subsequently, the streamlines continue to be 
pressed to the equator. Of course, we cannot talk about the behavior of the solutions far 
from the shock front, where our solution can yield a qualitatively incorrect result. However, 
we see that there is one tendency at small distances from the shock front: the bending of 
the streamlines toward the equator behind the shock. 

Translated by V. Astakhov 



Fig. 1. — A schematic image of the interaction between the pulsar wind and the interstellar 
medium. Because of the increase in the energy flux density in the wind toward the equator, 
the shock (i) takes the form of a torus. The streamlines bend at the shock front and are 
pressed to the equator. The region between the shock and the contact discontinuity (2) is 
fllled with relativistic particles emitting synchrotron radiation. The dotted line marks the 
r <^ Rn region, where the plasma flow may be considered to be quasi-steady 
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Fig. 2. — The distribution of the volume emissivity of synchrotron radiation in the plerion 
produced by the wind from the pulsar PSR 0531 + 21. The cross section of the shock is 
similar in shape to two contacting circumferences. The shock location at the equator is 
taken to be 14", in accordance with the observed location of the torus inner ring. 
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Fig. 3. — The flow geometry near the shock front. 
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Fig. 4. — A family of solutions to Eq. (A. 30) Solid line is the solution that corresponds to 
the absence of additional perturbations in the postshock flow. 
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Fig. 5. — The distribution of streamlines for a flow with a — 0.25. The shock location for 
isotropic {dashed line) and anisotropic {solid line) flows. The bending of the streamlines 
at the shock and the further bending of the flow toward the equator result in the flow 
compression near the equatorial plane. 



